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Abstract. A widespread, intuitive and computationally inexpensive method to 
|i analyze light guidance through waveguide bends is by introducing an equivalent 

straight waveguide with refractive index profile modified to account for actual 
waveguide bend. Here we revise the commonly used equivalent-index formula, 
• i-h . ending up with its simple extension that enables rigorous treatment of one- and 

two-dimensionally confined, uniformly bent waveguides, including tightly coiled 
i c"| microstructure fibers, curved ridge waveguides and ring microresonators. We also 

, Ph . show that such technique is applicable only to waveguides composed of isotropic or 

uniaxially anisotropic materials, with anisotropy axis directed perpendicular to the 
curvature plane. 

> 

To understand and to predict light wave behavior at waveguide bends was of 
major importance and interest to integrated- and fibre-optics community from around 
1970s and onwards. Today this interest is stimulated largely by two developments: 
(i) increasing the packaging density of integrated-optic circuits, backed by minimizing 
integrated waveguide bend radii while keeping bend losses at a tolerable level; and (ii) 



the advent of photonic crystal fibers possessing quite complicated, high index-contrast 
dielectric profiles as compared to step- or graded-index fibers for which early theoretical 
methods to treat bend losses were developed. 

Numerical methods to simulate light propagation through waveguide bends, such as 
^ beam propagation method pQ, the method of lines PIE], as well as more general-purpose 

finite-element and finite-difference techniques, can be very demanding computationally: 
in one recent example [4], a 64-processors cluster was used in full-vector finite-element 
modelling. No wonder one leans to analytic techniques to reduce computation burden 
whenever possible; one such technique, perhaps the simplest and most widely used 
in modelling of microstructure fibre bends today El 13 El, is equivalent-profile 
method (HlUDl- It reduces dimensionality of the problem by replacing actual waveguide 
bend by a straight piece with refractive index profile 

n 2 = e c (l + y/Rf ~ e c (l + 2y/R), (1) 

where ec is dielectric permittivity in the bent waveguide cross-section, as measured 
in Cartesian frame, R is the bend radius pointing from the x-directed curvature axis 
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to a (somewhat arbitrary) waveguide central plane, y the distance from that plane. 
Designed for elegant and inexpensive treatment of bends, though, formula © grounds 
on approximations whose validity may often be questioned, as those of weak-guidance 
regime and small curvature (y/R <C 1), and originally was introduced for quite a specific 
case of non-magnetic, step-index, low index-contrast fibre. In this Letter we overcome all 
these limitations by an alternative, first-principle derivation of expressions for modified 
dielectric permittivity e and magnetic permeability \i analogous to (0), but applicable 
to waveguides of arbitrary cross-section under the only assumption that bend radius is 
constant (while not necessarily large). 

Following Schouten jTT] and Post ^2] , we start from the generally covariant Maxwell 
equations 

29[a-E^] = — B Xu , d[ K B X u] = 0, (2) 
d u S) Xu = I) A + j\ <9 A £ A = p. (3) 

Here the square brackets denote alternation: •[>„] = ^{- Xu — va)> E\ and B\ v = —B v \ 
are covariant electric vector and magnetic bivector fields coinciding with electric 
and magnetic three-vectors E and B in Cartesian frame, $) Xl/ = —S) uX and D x 
are contravariant magnetic bivector density and electric vector density of weight +1 
corresponding to H and D, j A and p are electric current and charge densities, A, v — 
1 ... 3. With such transformation characteristics assigned to the electromagnetic field 
quantities, Eqs. (J2J), © are known to be form- invariant [H1II2], i>e., they do not change 
their form under arbitrary revertible coordinate transformations. It might be practical 
to convert Eqs. (J2J), © with use of dual equivalents <B K = \<£ kXv B Xv and H K = \I kXv $) Xu ', 
where <£ kXu and t K \ v are pseudo (hence tildes) permutation fields equal to Levi-Civita 
symbol in any coordinate system, to the form directly reminiscent of Maxwell equations 
in Cartesian frame: 

~^ Xv d x E v = -p KX h x , K p KX H x = 0, (4) 
fr^dxH,, = e KX E x + f, d K e KX E x = p, (5) 

while the constitutive relations are implied being 

<B A = /i Xu H v , £ A = e Xv E v . (6) 
Here e Xv and p Xv are tensor densities of weight +1: 

e Xu = |A|- 1 j£j£e AV , p Xv = |A|- 1 J£J> AV , (7) 

as stipulated by transformation characteristics of the fields H x , Q3 A , E x and D x . We 
denote by J x , = d x >x x the Jacobian transformation matrix, and A = det J x , its 
determinant. In this formulation, geometry enters Maxwell equations (j£jL (|5]l through 
material fields (0) exclusively; since the form of ©is precisely as if they were written 
in Cartesian components, a multitude of analytic and numeric methods developed for 
rectangular Cartesian frame apply. Furthermore, whenever e Xv and p Xv happen to 
be independent of one (or more) of the coordinates in a chosen coordinate system, 
that coordinate (or those coordinates) can be separated in the usual manner, like z 
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v = g 

owing to transformation behavior of the fundamental (metric) tensor g Xu 



coordinate in the case of straight homogeneous waveguide. Now we show that in 
cylindrical coordinates {x,r, <f)}, such ruling out of angular coordinate is the case 
for, in particular, isotropic guide of arbitrary cross-section, bent along 0, and specify 
Eqs. (J7J) for that case. 

Since isotropic media are defined as those possessing scalar (though, in general, 
position-dependent) permittivity €q and permeability /x c , as referred to Cartesian 
system, transformation rules (J7J) reduce to 

A" ,. — A., A'', ..A;' ., — ~.,A;<,. ^g^ 

T x T u n x ' u ' ■ 

its determinant g = det g Xv = det 2 Jy insofar as det g x ' v ' = 1 in Cartesian frame. 
Equations (JBJ) resemble 'effective' permittivity and permeability introduced in albeit 
we would refrain from using the word 'effective', on the ground that (jHJ) are nothing but 
transformation rules for the permittivity and permeability of isotropic media. With 
transformation from orthogonal Cartesian to cylindrical coordinates given by 

x — x, y = rcos<j), z — rsin<p. (9) 

(x defines bending axis), one can find the Cartesian-to-cylindrical transformation matrix 
for contravariant components 
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g-\ = r. After substituting these expressions in Eqs. (jHJ), separating the variable in 
Maxwell equations (HJ), © and differentiating the exp(i/30) multipliers with respect to 
0, normalizing the propagation constant /3 by R in the usual manner and introducing 
shifted coordinate y = r — R, one obtains 

e xx = e yy =e c (l + y/R), e zz = e c (l + y/R)' 1 (12) 

(we do not distinguish between co- and contravariant indices once particular coordinate 
system is chosen), and similarly for permeability components: 

fjvx = /% = vc{i + y/R), t^zz = nc{i + y/ R y 1 - ( 13 ) 

Note that expressions for e xx and e yy , multiplied by their magnetic counterparts, lead 
to refractive tensor components n\ x = e xx fi xx and n 2 yy = e yy fi yy precisely in line with 
standard formula © in the case of non-magnetic media. The e zz and [i zz components 
manifestly differ from the rest, however (see figure [IJ; this difference becomes especially 
pronounced when departing from weakly guiding approximation, as soon as non- 
negligible ^-components of electric and magnetic fields would 'probe' the e zz (x,y) and 
fi zz (x,y) profiles then. 
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A benchmark example to illustrate the validity of formulae (|12j) . (|13j) is a problem 
of light propagation through a homogeneous slab waveguide bend that permits precise 
analytic treatment — see e.g. ^3] for an overview. In our approach, covariant Maxwell 
curl equations, <£ KXu d\E u = —fi KX H\ and (t KXu d\H u = e KX E\, reduce to four independent 
first-order scalar equations then, grouped in two pairs: one for the E x , components 
(TE mode): 

d r E x = ik^H^, (14) 
r 

d r % = ik (e c r - 7 ^-\ E x , (15) 

where k is the free-space wavenumber, n e g = (3/(kR) the dimensionless mode index; and 
another pair for E^, H x (TM mode). To derive second-order equation for E x (which can 
be regarded as an eigenproblem in n off ), we differentiate (fT^j) with respect to r, use (fT3j) 
for drHcf, and ikH^ = -^d r E x to exclude H^. This leads, in the regions of constant /zc, 
to 

d 2 E x ldE x )2 / n 2 cS R 2 \ „ . . 

a Bessel equation used customarily in exact analysis of bent slab waveguides [HJ IT3] . 
while an attempt to get similar equation with use of approximate equivalent- index 
formula (0) fails. The results of full-vector two-dimensional finite-difference frequency- 
domain (FDFD) modelling of confined modes in a high-contrast step-index optical fibre 
(figure 2) demonstrate systematic discrepancy between approximate model (JTJ) and 
formulae (jl2j) . (fTHj) ; a discrepancy which can not be tolerated in realistic simulations of 
high index contrast ridge waveguides and holey fibers, ring microresonators, or sharp 
fibre bends such as those due to non-perfect alignment of fibre splices. 

A closer inspection of (fTT)|) brings to conclusion that the only type of initial 
anisotropy in ec and fic which still permits ruling out the <ft dependence in e Xv and fi Xu 
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Figure 2. The difference in fundamental mode indices at various bend radii R of a 
step-index fibre of radius o and dielectric index n core = f.45 in the air background, 
calculated numerically with and l|12Jl , . 



transformed according to (jZJ) is that of uniaxial crystal with anisotropy axis pointing in 
the x direction (i.e., along the bend axis): 

/ e b \ 

ec= e a . (17) 

V o o t a ) 

In that case, one gets instead of formulae 1)120: 

£^ = 66(1 + 3//^), (18) 
e yy = e a (l + y/R), (19) 
e z , = e (l + y/i?)" 1 , (20) 

and likewise for modified Unfortunately, however, the case of anisotropy given by (|T7|) 
is not one encountered when modelling, e.g., photonic bandgap fibre filled with liquid 
crystal; for that and other examples of curved waveguides possessing anisotropy of some 
general kind, rigorous equivalent-profile formulation is not applicable and should be 
substantiated by more brute-force numeric or less stringent analytic techniques. 

In summary, equivalent-profile expressions (J12)L (|T3*j) to treat waveguide bends 
have been presented. Unlike conventional formula (fl)). they hold for arbitrarily tight 
bends, arbitrary high-contrast material profiles, and permit the inclusion of magnetic 
parts. This offers three-dimensional modelling accuracy with those existing full-vectorial 
two-dimensional finite-difference or finite-element solvers which permit 3x3 diagonal 
matrices for e and \x be put in. Eliminating the restraints on bend radius and refractive 
index contrast enables one to simulate also the whispering-gallery modes in spherical, 
toroidal and other types of optical microcavities. The principal limitation of the 
equivalent-profile technique is that many cases of non-trivial anisotropy can not be 
treated rigorously in the manner above. 
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